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Abstract 

A representation for integral kernels of Delsarte transmutation operators is obtained in the 
form of a functional series with the exact formulae for the terms of the series. It is based on 
the application of hyperbolic pseudoanalytic function theory and recent results on mapping 
properties of the transmutation operators. 

The kernel K\ of the transmutation operator relating A — — j-^ + qi (x) and B = — results 
to be one of the complex components of a bicomplex-valued hyperbolic pseudoanalytic function 
satisfying a Vekua-type hyperbolic equation of a special form. The other component of the 
pseudoanalytic function is the kernel of the transmutation operator relating C = — -j^ + 12 (x) 
and B where 92 is obtained from q± by a Darboux transformation. We prove the expansion 
theorem and a Runge-type theorem for this special hyperbolic Vekua equation and using 
several known results from hyperbolic pseudoanalytic function theory together with the recently 
discovered mapping properties of the transmutation operators obtain the new representation for 
their kernels. Several examples are given. Moreover, based on the presented results an approach 
for numerical computation of the transmutation kernels is proposed. 



1. Introduction 

The notion of a transmutation operator relating two linear differential operators was 



introduced in 1938 by J. Delsarte 12 and nowadays represents a widely used tool in the 
theory of linear differential equations (see, e.g., [l], [7j, [32| , 34 , 42 , |44| ) . Very often in 



literature the transmutation operators are called the transformation operators. Here we keep 
the original term introduced by Delsarte and Lions [13] . It is well known that under certain 
regularity conditions the transmutation operator transmuting the operator A = —^2 + q(x) 
into B = —^2 is a Volterra integral operator with good properties. Its integral kernel can be 
obtained as a solution of a certain Goursat problem for the Klein-Gordon equation with a 
variable coefficient. There exist very few examples of the transmutation kernels available in a 
closed form (see |26| ). 

In the present work we obtain a representation for the kernels of the transmutation operators 
for regular Sturm-Liouville operators (with complex valued coefficients) in the form of a 
functional series with the exact formulae for the terms of the series. The result is based on 
several new observations. We use our recent result on the construction of the kernel of the 
transmutation operator corresponding to a Darboux associated Schrodinger operator 26 to 
find out that a bicomplex-valued function whose one complex component is the transmutation 
kernel Ki(x,t) (for a Schrodinger operator — qi(x), q x e C[—b,b]) and the other complex 
component is K%(x,t) (for a Schrodinger operator — Q2(x), with q2 obtained from q\ by a 
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Darboux transformation) is a solution of a certain hyperbolic Vekua equation of a special form 
(for the theory of elliptic Vekua equations see 45 as well as [2] and |23| ). 
In spite of a recent progress reported in [20) . |23| . 25 , |29| the theory of hyperbolic Vekua 



equations is considerably less developed than the theory of classical (elliptic) Vekua equations. 
For example, as was shown in |25| (see also [23] ) the construction of so-called formal powers 
(solutions of the Vekua equation generalizing the usual complex powers [z — z ) n ) can be 
performed using the definitions completely analogous to those introduced by L. Bers [2]. 
Nevertheless no expansion theorem nor a result on the completeness of the obtained system 
(a Runge-type theorem) of hyperbolic formal powers was available. In this work we apply the 
results from [25| for constructing the formal powers for the hyperbolic Vekua equation arising 
in relation with the kernels K\ and Ki as well as the observation from [fjj establishing that they 
are the result of the transmutation of the usual powers of the hyperbolic variable. We obtain an 
expansion theorem and prove a completeness result. Moreover, we obtain explicit formulae for 
the expansion coefficients which leads to the functional series representation for the kernel K\ 
(and also for K%). We give several examples of application of this result including a numerical 
computation. Finally, we propose an alternative method for approximate construction of the 
transmutation kernel based on the same hyperbolic formal powers but instead of the expansion 
theorem using the obtained completeness result. 

The paper is structured as follows. In Section 2 we introduce the hyperbolic and bicomplex 
numbers, the Vekua equation one of the solutions of which is the bicomplex function K\ — ]K%. 
We present several properties of that Vekua equation, construct an infinite system of its 
solutions called formal powers, explain its relation to generalized wave polynomials and 
to spectral parameter power series (SPPS) representation for solutions of Sturm-Liouville 
equations. We consider a Goursat problem for the hyperbolic Cauchy-Riemann system and 
obtain certain related results important for what follows. In Section 3 we recall some long known 
facts about the transmutation operators and give some recent and new results concerning the 
relation between K\ and Ki , the commutation and the mapping properties of the corresponding 
transmutation operators. In Section 4 we introduce an operator transforming the boundary 
data of the Goursat problem for the wave equation into the boundary data of the Goursat 
problem for the operator □ — q{x). Using its properties we obtain a result about a functional 
series representation for solutions of the equation (□ — q) u = in terms of generalized wave 
polynomials. In Section 5 we introduce a transmutation operator relating the hyperbolic 
Cauchy-Riemann operator with the interesting for us hyperbolic Vekua operator. The main 
observation here is that this transmutation operator possessing important boundedness 
properties maps usual powers of the hyperbolic variable z into the formal powers as well 
as the usual derivatives with respect to z into the generalized derivatives in the sense of 
L. Bers. Here we mainly use the ideas from [5] where similar observations were made in 
the elliptic setting. Section 6 is dedicated to the expansion and the Runge-type theorems 
corresponding to the special Vekua equation under consideration. We show that a generalized 
Taylor formula for its solutions is valid under certain conditions on their Goursat data. The 
formula involves the formal powers constructed and studied in preceding sections. The obtained 
Runge-type theorem establishes that even when the conditions required for the validity of the 
Taylor expansion are not fulfilled any solution of the studied hyperbolic Vekua equation can be 
approximated arbitrarily closely by a formal polynomial. In Section 7 we show that K% and K 2 
are necessarily the conjugate components of a solution of the Vekua equation. Using this we 
give explicit formulae expressing Ki in terms of K\ and vice versa. Section 8 contains the main 
result of this work - the representation of the transmutation kernel in terms of a generalized 
Taylor series in formal powers. For the corresponding expansion coefficients we obtain both a 
recurrent and a direct exact formulae. Several examples are given illustrating their application. 
Finally, in Section 9 we show how the obtained results can be used for numerical computation 
of the transmutation kernels. Here together with the approach based on the generalized Taylor 
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expansion theorem we propose a method based on the Runge-typc theorem. When sufficiently 
many derivatives of the potential taken in the origin are not available this second method 
is more convenient. A developed numerical example illustrates its excellent computational 
performance. 



2. Elements of hyperbolic pseudoanalytic function theory 

2.1. Hyperbolic and bicomplex numbers 

By j we denote the hyperbolic imaginary unit: j 2 — 1 and consider the algebra of hyperbolic 
numbers, also called duplex numbers (see, e.g., [43) ) 



D := {c = a + bj : j 2 = 1, a, b e R} = C1 R (0, 1). 
The algebra D is commutative and contains zero divisors. Additional information on the 



hyperbolic numbers can be found in 31 , 36 and 43 



We will consider the variable z — x + t j where x and t are real variables and the corresponding 
differential operators 

d z = - (d x +jd t ) and = - {d x - jd t ) ■ 

As in the case of complex numbers, we have d s z = which explains the choice of the minus 
sign in the definition of d g . 

Let B denote the algebra of bicomplex numbers which can be defined as follows 

B := {w = u + v j : it, v E C} 

and the complex imaginary unit i commutes with j. More on bicomplex numbers can be found, 
e.g., in [I], (I], [23] and [Si 



The conjugate of a bicomplex number w with respect to j we denote by w, i.e., w = u — vj. 
The corresponding conjugation operator is denoted by C, Cw = w. 

C-valued functions will be also called scalar. For the scalar components of w we introduce 
the notations 

lZ(w) =u = -(w + w) and I(w) = v = —{w — ID). 
Z zj 

In order not to overload the text by excessively many different notations we will use the 
notation 1Z and X also in the operational sense as projection operators projecting aieB onto 
the respective scalar components, TZ =j(7 + C) and 1= ^(J— C) where I is the identity 
operator. 

It is convenient to introduce the pair of idempotents P + = |(1 + j) and P~ = |(1 — j) 
satisfying (P ± ) 2 = P ± and P + P~ = P~P + = 0. Then for any w <E B there exist the unique 
numbers w + , w~ G C such that w = P + w + + P~w~ which are related with the components 
of w in the following way 

w ± =TZ{w)±2(w). (2.1) 

A nonzero element w £ B belongs to the set of zero divisors <r(B) iff w = P + w + or w = P~w~ . 
We will use the following norm (see [4]) in B, 

M = 5(KI C + Mc)> ( 2 - 2 ) 

where l-L is the usual norm in C. 
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2.2. Hyperbolic Vekua equation 
We will consider functions w : I 



and Vckua-type equations of the form 



d~w = 



bw 



(2.3) 



where a and b are bicomplex functions of the variable z G D. In 20 , 23 , 25 it was shown that 



many results from pseudoanalytic function theory [2] remain valid in the hyperbolic situation 
We will not give here definitions and properties corresponding to the general Vekua equation 



d- z W 



-W = 



(2.3) referring the reader to [23] . Instead, as we are interested in a very special form of the 
coefficients a and b we restrict ourselves to that particular case and obtain several results, such 
as the expansion theorem and a Runge-type theorem which are still unavailable in a general 
situation. Thus, let us consider the following hyperbolic Vekua equation 

I 

where / is a scalar, non-vanishing function in the domain of interest. Vekua equations of this 
form, i.e., with a = and b being a logarithmic derivative of a scalar function are called Vekua 
equations of the main type or main Vekua equations. They arise from the factorization of the 
stationary Schrodinger operator (in the elliptic case) 
Gordon operator (in the hyperbolic case) 



(2.4) 



21 , 23 and of the variable mass Klcin- 



23 . Moreover, for the purposes of the present 



paper it is sufficient to consider the case when / depends only on x and hence (2.4) can be 
written also in the form 

/'(*) 



d- z W 



It is equivalent to the system 



fd x 



2/(x) 



w = 0. 



d t v, 



1 



dx {fv) = d t u, 



f 

where u, v: M 2 — > C and W = u + Introducing 
rewrite this system also in the following form 

1 



f 



and 



the notation <p — uj f and ip 



(2.5) 

(2.6) 

(2.7) 

fv one can 

(2.8) 



In the case when all the involved functions were real valued this system was studied in |37 



38 and later on in 25 , 23 



Equation (2.4) admits a corresponding generating pair F — f and G — j/ /. We recall that 



a generating pair corresponding to a Vekua equation is a pair of its solutions F and G, 
independent in the sense that an arbitrary function w, in this case B-valued, can be represented 
in the form 



uj = ipF + tpG 



(2.9) 



where ip and 4> are scalar functions. That is, F and G should satisfy the condition I(FG) ^ 
everywhere in the domain of interest. In relation with the generating pair (F,G) = (/,j//) 
observe that the scalar functions ip and ip satisfy (2.8) iff the B-valued function uj — ipf + ipj/ f 



is a solution of ( 2.4 ) 



The knowledge of a generating pair allows one to define the (F, G)-derivative in the sense of 
L. Bers. For a function to written in the form (2.9 1 its (F, G)-derivative has the form 

(d z <p)F+(d z ^)G. 



Z (F,G) 

dz 



CONSTRUCTION OF TRANSMUTATION OPERATORS 



5 



Vekua equation 23 



Whenever W satisfies ( |2.4| its (F, G)-derivative (with (F,G) = (/, j//) ) is a solution of the 

d z f. 



Notice that for / = f(x) this equation can be written as follows 

d- z w + i^\w = 



(2.10) 



or 



2f(x) 

(f-Hx))' 



w 



0. 



(2.11) 



Thus, the coefficient in the equation has the same structure as in ( |2.5[ ) but the role of / is 
played by 1//. Hence a generating pair for (2.11) can be chosen as {F\,G\) = (1//, /j). The 

f I 00 o ° 

second derivative of a solution of (2.5 1 can be defined in the form W = w where w = W = 
fiFi + i/jiGi = (pi/ f + ipxfi, <pi and ipi are scalar functions, and w = (d z <pi) If + (d z ipi) /j- 
Obviously, the next generating pair can be chosen again as [F2, G2) = (/, j//) and in this way 
the derivative of any order in the sense of Bers can be defined. The derivative of n-th order we 
denote as W^. In the theory of pseudoanalytic functions such a sequence of generating pairs 
is known as a periodic generating sequence (of the period two). We have (F n ,G n ) = (/J//) 
when n is even, and (F n , G n ) = (1//, /j) when n is odd. 



For what follows it is important to observe that if W is a solution of (2.5 1, its n-th derivative 
in the sense of Bers (if exists) has an especially simple form 



: j™9 t "W. 



(2.12) 



Indeed, for W = <pf + ipj/f we have W = {d z ip) / + (d z tp) j// where ip and ip are solutions of 
EM. Consider 



d z f = ^ (d x (p+jdt<p) = ~ 



1 



and similarly, 



1 



\ (fd t v+jd t i>) - f -d t 



ipf + ip 



= -d t W. 
2 



Thus, W = jdtW. As this reasoning does no t chan ge if / is substituted everywhere by 1// we 



obtain that W — jdtW = j 2 d^W and hence ( 2.12 ) is valid both for odd and for even values of 



Below we extend some of the results on the relationship between hyperbolic pseudoanalytic 
functions and solutions of the Klein-Gordon equations to the bicomplex case. Since all the 
proofs are essentially the same, we refer the reader to |23| Chapter 13], |25| for details. The 
operator <9 2 applied to a scalar function ip can be regarded as a kind of gradient. If dzf = 
where $ is a B-valued function defined on a simply connected domain with $1 = 7?.($) and 
$2 = such that 

5 t $i+a x $ 2 -0, V(i,t)6n, 

then we can construct tp up to an arbitrary complex constant c. Indeed, let T C O be a rectifiable 
curve leading from (xo,to) to (x,t), then the integral 



A h [$]{x,t) :=2 



$1 dx - $2 dt 



(2.13) 



G 
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is path- independent, and all C-valued solutions tp of the equation dgip = $ in f2 have the form 
cp(x, t) — Ah[<&(x, t)] + c where c is an arbitrary complex constant. In other words the operator 
Afi denotes the well known operation for reconstructing the potential function from its gradient. 
In Section[7|we need a particular case of (2.13) when V is the path consisting from 2 segments, 
first going from (a;o,io) to (x,to) and second going to (x,t). Assuming that this path belongs 



to the domain of interest fl, in such case formula (2.13) reads as follows 



A h [*](x,t) = 2 



<f>i(r],t )dri 



to 



(2.14) 



The following theorems explain relationship between the main Vekua equation and variable 
mass Klein-Gordon equations. We use the notation □ := d 2 — 9 t 2 . 



Theorem 2.1 [25 . Let f = f(x) be a non-vanishing scalar function and W be a solution 

(2.15) 



of (2.5). Then u = 1Z(W) is a solution of the equation 

(□ — qi(x))u = 0, where qi 
and v — Z(W) is a solution of the equation 



r 



(D — q 2 (x))v = 0, where q 2 = 



— o 

f- 1 v/ 



(2.16) 



Theorem 2.2 25 . Let u be a scalar solution of the Klein-Gordon equation (2.15) in a 



simply connected domain Q. Then a scalar solution v of the associated Klein-Gordon equation 



(2.16) such that u + ju is a solution of (2.5) in $1, can be constructed according to the formula 



v = -r i A h [ifd- z (r i u)]. 

It is unique up to an additive term cf~ l where c is an arbitrary complex constant. 



Vice versa, given a solution v of (2.16), the corresponding solution u of (2.15) such that 



u + ji> is a solution of (2.5), has the form 

u=-fA h [if- 2 d,(fv)]., 

up to an additive term cf. 



2.3. Integration and formal powers 



Let (F, G) be a generating pair. Its adjoint generating pair (F, G)* = (F* , G*) is defined by 
the formulae 

2F G , 2G 



FG-FG FG- FG 

For our special case of generating pairs we have (/J//)* = (/j, 1/7) and (1//, /j)* = (j//, /). 
The (F, G)-integral is defined as 



Wd 



1 



z = ~[F(z 1 )K 



G*Wdz +G( Zl )K 



F*Wdz 



where T is a rectifiable curve leading from zq to z\. 

If W = tpF + ipG is a bicomplex (F, G)-pseudoanalytic function where ip and tp are scalar 
functions then 

Wd {FtG) z = W{z) - <p(z )F(z) - tp(z )G(z), 
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and this integral is path-independent and represents the (F, G)-antiderivative of W. 

Definition 2.3 [2]. The formal power Zm\a, zq;z) with center at Zq, coefficient a and 
exponent is defined as the linear combination of the generators F m , G m with scalar constant 
coefficients A, /i chosen so that XF m (zo) + (j,G m (za) = a. The formal powers with exponents 
n = 0,1,2,... are defined by the recursion formula 



z£ +1 \a 1 z a ;z) = (n + l) 



7 {n) 
J m+ 



i(a,z -X) d (F m ,G m )C- 



This definition implies the following properties. 

(1) Z<S\a,Zo; z) is an (F m , G m )-bicomplex pseudoanalytic function of z, i.e., it satisfies the 
hyperbolic Vekua equation determined by (F m ,G m ) (see |23| ). 

(2) If a' and a" are scalar constants, then 

+ jo", ^o; z) = o'zW(l, z ; z) + o"zW(j, z ; z). 

(3) The formal powers satisfy the differential relations 



d(F m ,G m )Zm\a 1 z a ;z) 



dz 



nZ { ™ + P(a,z Q ;z). 



(4) The asymptotic formulae 



z m\ a > z o; z) ~ a(z - z ) n , z 



zo 



hold. 

In the following we denote Zq (a, z ; z) by Z^ n \a, z Q ; z). 

Since we are primarily interested in the hyperbolic Vekua equation (2.51, observe that as was 
explained above, a corresponding generating sequence is periodic having the form (F n , G n ) = 
(fjj/f) when n is even, and (F n ,G n ) — (1//, /j) when n is odd. Using the construction from 
[2] given there for an analogous elliptic situation (see also |23[ Sect. 4.2]) we obtain that 
the corresponding formal powers admit the following elegant representation. We consider the 
formal powers with the center in the origin and for simplicity assume that /(0) = 1. Define the 
following systems of functions 



X^{x)=X^(x) = 1, 



X {n \x)=n 
X {n \x)=n 



ds, 



X 



(n-l) 



(-1)" 



Then the formal powers corresponding to (2.5| have the following form. For a = a' 
have 

j 



Z^(a,0;z) = f(x)H(,zW{a,0;z)) + -J-rl(,Z^(a, 0; z)), 

f(x) 



(2.17) 
(2.18) 

(2.19) 



■ja" we 



(2.20) 



where 



fe=0 



Z^(a,0;z) =a'J2 "U M) jV +3 a "J2 (?p (n_fc) j fc ** for an odd n (2.21 



fe=0 



and 



*ZW(a,0;z) =a'J2 (^X^-^ft* ff\ X ^-^j k t k for an even n. (2.22 

k=0 ^ ' fe=0 ^ ' 

For any a € B and n € N = N U {0} the formal power (|2.20[) is a solution of \2.ty. 
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Remark 2.4. Formulae ( 2.20 1— ( 2.22) clearly generalize the binomial representation for the 
analytic powers az n . If one chooses / = 1 then Z^ n \a, 0; z) = az n . 

The formal powers can be written with the use of generalized wave polynomials 



and (2.191 by 



18 . Consider two families of functions {<Pfc}^ and {ipk}kLo constructed from (2.17), (2.18) 



' f{x)X {k \x) 1 fcodd, 
f(x)X { - k) (x), k even, 



(2.23) 



and 



i>k{x) 



/(f) 



k odd, 
k even. 



(2.24) 



The following functions 

Uo=<Po(x), «2m-l(M) 

and 



fe=0 



odd fe=l 



W O =-0o(2;), «2m-l(^,*)= J! ( uj^m-kix)^, V 2m (x,t) = ^ 



even fc— 



odd k=l 



V>m-k(x)t k , 

(2.25) 



(2.26) 

are called generalized wave polynomials (the wave polynomials are introduced below, in Remark 
2.9). It is easy to see that 



Z (o) (a,0;z) = a'u (x,t) + ja"v (x, t), (2.27) 
2 (n) (a,0;z) = a'u 2n -i(x,t) + a"u 2n (x,t) + j(a'v 2n (x,t) + a"v 2n -i(x,t)), n>l. (2.28) 

The following parity relations hold for generalized wave polynomials. 

u (x, -t) = u (x,t), u 2n -i(x,-t) = u 2n -i(x,t), u 2n (x, -t) = -u 2n (x,t). (2.29) 

It is worth to mention the following result obtained in [22] (for additional details and simpler 
proof see 23 and 24 ) which establishes the relation of the systems of functions {^fclfcLo anc ^ 
{V'felfcLo t° the Sturm-Liouville equation. 

Theorem 2.5 [22], Let q be a continuous complex valued function of an independent real 
variable x <E [a, b] and A be an arbitrary complex number. Suppose there exists a solution f of 
the equation 

f" -qf = (2.30) 

on (a, b) such that f g C 2 (a, b) D C 1 [a, b] and f(x) ^ for any x € [a, b]. Then the general 
solution g g C 2 (a, b) D C 1 [a, i>] of the equation 



on (a, 6) has the form 



.9" - 19 = A.9 



5 = cigi + c 2 g 2 
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where c\ and c 2 are arbitrary complex constants, 



(2fc + i; 



(2.31) 



and both series converge uniformly on [a, b] together with the series of the first derivatives 
which have the form 



g'l = f' + J2 ToTTTT ( T^k + 2fc^ 2fe _! 



1 (2^)! V / 



and 



00 \k 
■92 = E ^ 



fc=0 



(2fe + l)! V/ 



<^ 2fc+ i + (2k + 1) V 2 fe ■ (2.32) 



The series of the second derivatives converge uniformly on any segment 
[01,61] C (a, 6). 



Representations (2.31 1 and ( 2.32[ ), also known as the SPPS method, present an efficient and 
highly competitive technique for solving a variety of spectral and scattering problems related 
to Sturm-Liouville equations, see [§], |14| , 17 , 19 , 24 and 27 and references therein. 



2.4. The Goursat problem for the hyperbolic Cauchy-Riemann system and a Runge type 
theorem 

For the construction of the kernels of the transmutation operators we are interested in solving 
corresponding Goursat problems. Our results concerning the Goursat problems for equation 



(2.5 1 arising as a natural step in that construction are based on the following simplest case. 
Consider the hyperbolic Cauchy-Riemann system 

dziu = 



u t = v x 



(2.33) 



studied in a considerable number of publications (see, e.g., 31 , 36 , 47 ). It is easy to see 
that its general solution can be written as follows 



w(x, t) 



x - 1 



-j * 



x + 1 



x — t 



x + t 



x - t 



(2.34) 



where $ and VE" are arbitrary continuously differentiable scalar functions. 

Considering the values of w on the characteristics x — t and x — — t we see that 

w(x,x) = P + <$>(x) + P~4'(0) and w(x, -x) = P + $ (0) + (x) . 



(2.35) 



Thus, we arrive at the following correct statement of the Goursat problem. Find a solution 



w of (2.33) in a closed square R with the vertices (±26,0) and (0,±26), 6 > satisfying the 



conditions (2.351 where $ and ^ are arbitrary scalar functions such that e C 1 [— 6, b]. 



Proposition 2.6. The unique solution of the Goursat problem (2.331, (2.35) has the form 
(p04l 



Proof. The fact that w is a solution of (2.33) follows from (2.34). The conditions (2.35) 
are obviously satisfied. To prove the uniqueness one considers the Goursat problem with 



10 



VLADISLAV V. KRAVCHENKO AND SERGII M. TORBA 



the homogeneous conditions and finds out that $ and \1/ from the representation ( 2.34[ ) of 
the general solution of ( |2.33| satisfy the equations P + $(x) + P~*(0) = and P+<1>(0) + 
P~'$>(x) = 0, x e [—6, b] from which both $ and *Sf are necessarily trivial. □ 

This proposition is not entirely new. As a kind of analogue of the Cauchy integral formula 



for analytic functions it was discussed, e.g., in 36 in the case w : D — > D. For our purposes it 



is more convenient to consider this simple representation of any hyperbolic analytic function 
via its values on the characteristics from the point of view of the Goursat problem. 



The following result is a direct corollary of Proposition 2.6 



Proposition 2.7. Let $ and 4" be scalar functions defined on [—b, b] and admitting there 
the uniformly convergent power series expansions 

oo oo 

$(x) = ^a n x n and 3> (x) = ^ j3 n x n . 



71=0 



n=0 



Then the unique solution of the Goursat problem (2.331, (2.35) has the form 



w 



n=0 



where 



(P+a n + P~l3 n ) 



(2.36) 



(2.37) 



Proof. Due to Proposition |2.6| we have 



,t)=P+J2 a n 



ri=0 



X + t 



n=0 



X — t 



n=0 ^ ' 



Notice that P ± z = P ± (x± t) and hence P ± z n = P ± (x ± t) n . Thus, 



E + ^P~A =E4 (P + -n + P-Pn) 



ri=0 



n=0 



□ 



Remark 2.8. Under the conditions of Proposition |2.7| we have that the unique solution 



Taylor series (2.361 where 



w of the Goursat problem ((2733 ), (2.351 can be represented as a uniformly convergent in R 

_ iy[™](0) 



(2.38) 



and due to (2.34), 



Indeed, let us prove this formula for the coefficients. Due to (|2.12| we have that w^ n ' — j™9™u> 

ra> =^ (p+p<£ {n) (^) + p-f (-i)" ( 1 ' 

Now we notice that P+j n = P+ and p—f (-1)" = P~ for any n = 0, 1 Thus, w^(0) 

£(P+a n + p-p n )=n\a n . 
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Remark 2.9. It is possible to give another representation for the solution of the Goursat 
problem, namely, a representation as a series in wave polynomials, cf. 18 Proposition 1]. 
Recall that the following polynomials 

p (x,t) = l, P2m -i(x,t)=K({x+jt) m ), p 2m (x,t) =T((x+it) m ), m> 1, 

are called wave polynomials. They may be also written as 



p (x,t) = l, P2m-l(x,t)= ^ 



even k—0 



X m ~ k t\ P2m (x,t)= 

odd k=l 



x m - k t k , (2.39) 



cf. with the definition of generalized wave polynomials (2.25), (2.26). Since z™ = P2n-i( x ,t) + 
JP2n(x,t) for n > 1, we obtain from (2.37) that 

w(x, t) = (P+a Q + p-/3 )p (x, t) 

™P+a n + P-/3 n ^ P+a n -P-(3 n 

+ 2^ P2n-l{X,t) + 2^ ^ P 2n(X,t). 



n=l 



n=l 



2" 



Obviously, an arbitrary continuously differentiable solution of (2.33) cannot be represented 



globally as a uniformly convergent power series of the form (2.36). Nevertheless the following 



analogue of the Runge approximation theorem from classical complex analysis is a corollary of 
the preceding proposition. 



Proposition 2.10. Let we C 1 (R) be a solution of (2.33) in R. Then there exists a 

N _ 

sequence of polynomials P^ — a n z " uniformly convergent to w in R. 

n=0 



Proof. We need to prove that for any e > there exist such a number N and coefficients a n , 
n = 0, 1, . . . ,N that \w(x, t) — Pjsr(x, f)| < e for any point (x, t) £ R. Let the scalar functions 
$ and ^ defined on [—b, b] be such that (2.35 1 hold. As was shown above such scalar functions 
always exist and under the condition of the proposition they are obviously from C 1 [— b, b]. We 
choose £ > and such £1,2 > that e — (e% + £ 2 )/2- According to the Weierstrass theorem there 
exists such a number N and such polynomials p\ and p 2 of order not greater than a certain TV 
that |$ (x) — pi(x)\ c < E\ and \^ (x) — p 2 (x)\ c < e 2 , —b<x< b. Due to Proposition 



unique solution w of the Goursat problem for equation (2.33) with the boundary conditions 



2.7 



the 



w{x, x) = P + px (x) + P p 2 (0) , w(x, -x) = P + pi (0) +P p 2 (x) 

has the form w = Pn where P N — Yl n =o §s* (P +a n + P fin) % n with a n and f3 n being the 
coefficients of the polynomials p\ and P2 respectively. 
Consider 



\w(x,t) — w(x, t)\ — \w(x,t) — Pn(x, t)\ 

'x + t 



2 

x + 1 



x — t 



-P+(x,i) 



-P+P+(x 7 t)-P-P N (x,t) 
x - t 



-P^(x,t) 



2 

where we have used tfHty and ( flE2j ). Notice that P£(x,t) = J2n=o a ™ (^)™ =V\ (^) and 

PN(x,t) = En= Pn(^) n = P2 (W TllUS, \ W (x, t) - P N ( X , t)\ < \{s X + £ 2 ) = E. ~ □ 
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3. Transmutation operators 
We give a definition of a transmutation operator from |27| which is a modification of the 



definition given by Levitan 32 adapted to the purposes of the present work. Let £ be a linear 
topological space and E\ its linear subspace (not necessarily closed). Let A and B be linear 
operators: E\ — > E. 



Definition 3.1. A linear invertible operator T defined on the whole E such that E\ is 
invariant under the action of T is called a transmutation operator for the pair of operators A 
and B if it fulfills the following two conditions. 

(i) Both the operator T and its inverse T _1 are continuous in E; 

(ii) The following operator equality is valid 

AT = TB (3.1) 

or which is the same 

A = TBT' 1 . 



Our main interest concerns the situation when A = — -jzz +q{ x ): B = — 4^2, and q is a 
continuous complex-valued function. Hence for our purposes it will be sufficient to consider 
the functional space E = C[a, b] with the topology of uniform convergence and its subspace 
Ei consisting of functions from C 2 [a, b]. One of the possibilities to introduce a transmutation 
operator on E was considered by Lions [33] and later on in other references (see, e.g., |34| ), 
and consists in constructing a Volterra integral operator corresponding to a midpoint of the 
segment of interest. As we begin with this transmutation operator it is convenient to consider a 
symmetric segment [—b, b] and hence the functional space E — C[—b, b]. It is worth mentioning 
that other well known ways to construct the transmutation operators (see, e.g., 



32 



44 



imply imposing initial conditions on the functions and consequently lead to transmutation 



operators satisfying (3.1 ) only on subclasses of E\. 

Thus, consider the space E = C[—b, b]. In [6] and |26| a parametrized family of transmutation 
operators for the defined above A and B was studied. Operators of this family can be realized 
in the form of the Volterra integral operator 



T h u(x) = u{x) 



K(x,t; h)u{t)dt 



(3.2) 



where K(x, t; h) = H(^±^ , h), h is a complex parameter, \t\ < \x\ < b and H is the unique 
solution of the Goursat problem 

d 2 U{u,v;h) 



du dv 
H( u ,(U) = J+2 



q{u + v)H(u, v; h), 



n(Q,v;h) = ~. 



(3.3) 



(3.4) 



If the potential q is continuously differentiable, the kernel K itself is a solution of the Goursat 
problem 

^-q(x))K(x,t;h) = ^ s K(x,t;h), (3.5) 



dx 2 



dt 2 ' 



K(x, x; h) 



h 
2 



q(s) ds, K(x, —x; h) 



h 
2' 



(3.6) 



If the potential q is n times continuously differentiable, the kernel ~K.(x,t;h) is n+1 times 
continuously differentiable with respect to both independent variables. 
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Remark 3.2. In the case h = the operator coincides with the transmutation operator 

it was established that in the case q € 



33 



44 



studied in [34j Chap. 1, Sect. 2]. In [32 

C 1 !— 6, 6] the Volterra-type integral operator (3.2) is a transmutation in the sense of Definition 



3.1 



on the space C [— b, b] if and only if the integral kernel K(a;, t) satisfies the Goursat problem 



(3.51, (3.6) 



The following proposition shows that it is sufficient to know the transmutation operator 
T/jj or its integral kernel ~K.h t for some particular parameter h\ to be able to construct 
transmutation operators or their integral kernels for arbitrary values of the parameter 
h. 



Proposition 3.3 [6], [27] . The operators T hl and T h2 are related by the expression 

h 2 - hi 
2 



T h .u = T 



hi 



u(x) 



u(t) dt 



valid for any u G C[— 6, b]. 

The integral kernels K(a;, t; hi) and K(x, f; /i 2 ) are related by the expression 



K(x, t- h 2 ) = + K(x, t-M)+ 



(K(x, s; hi) — K(x, — s; hij) ds. 



The following theorem states that the operators are indeed transmutations in the sense 
of Definition 13.11 



Theorem 3.4. Let q € C[—b, b]. Then the operator T h given by ( |3.2[ ) satishes the equality 

(-- 
\ dx 2 

for any u £ C 2 [—6, b] . 



2 +q(x) ) T h [u] = T h 



d 2 
dx 



(3.7) 



Remark 3.5. This theorem under additional assumptions on q was proved in [6| and [26] 



In 27 it was shown that (3.7) holds for any h whenever it holds for some particular value h\. 



Proof. It follows from Proposition 3.3 that it is sufficient to prove i3.7l for one particular 



h, see 27 Proof of Theorem 5.6] for details. We take h = and the operator T := Tq. Let 



K denotes the integral kernel of T. There exists a sequence of potentials {gVijneN C C 1 [— b, b] 
such that q n — > q, n — > oo uniformly on [—6,6]. Let K 9n and T qn be the integral kernels and 
the transmutation operators corresponding to the potentials q n . Then K 9ii — > K uniformly and 
T„ — > T in the operator norm, see, e.g., 34 Chap. 1, Sect. 2}. Take a function u £ C 2 [—b, b}. 
For the operators T gn the following equality holds 26 Theorem 6] 

d 2 



dx- 



+ q n (x) T qn [u 



dx' 



(u) 



which with the use of ( |3.2[ ) may be rewritten as 
dx 2 



K qn (x,t)u(t)dt = -u"(x) + q n (x)T qn [u](x) + T gn [u"](x). 



Let y n (x) := J_ K.q n (x, t)u(t) dt and z n 



§ x _ x K(x, t)u(t) dt and 



q n T qn [u] 



Then y n -> y 



z, n — > z := —u 



?TM+T[/], 



oo. Moreover, functions y n satisfy 
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the initial conditions y n (0) = y' n (0) = 0. Now (3.7) follows from the f act t hat the operator d 2 
with the domain {y 6 C 2 [— b, b] : y(0) = y'(0) = 0} is closed, see, e.g., |l6| . □ 

The following theorem summarizes the mapping properties of the operators T/j est ablish ing 
that there exists a value of the parameter h such that T/j maps x k to ipk defined by (2.23). 



Theorem 3.6 ]6 , 26 . Let q be a continuous complex valued function of an independent 
real variable x£ [—b,b] for which there exists a particular solution f of (2.30) such that 
f £ C 2 [~b, b], f ^ on [—b,b] and normalized as /(0) = 1. Denote h := f'(0) £ C. Suppose 



T/j is the operator defined by (3.2) and ipk, k £ No are functions defined by (2.23). Then 

T h x k = ip k (x) foranyk£N . 



(3.8) 



Moreover, T/j maps a solution v of an equation v" + uj 2 v = 0, where ui is a complex number, 
into a solution u of the equation u" — q(x)u + u> 2 u = with the following correspondence of 
the initial values 

u(0)=v(0), u'(0)=v'(0) + hv(0). 



Remark 3.7. The mapping property (3.8 1 of the transmutation operator allows one to see 



that the SPPS representations (2.31) from Theorem 2.5 are nothing but the images of Taylor 
expansions of the functions cosh \J Xx and ^ sinh VAx. Moreover, equality (3.8) is behind a 



new method for solving Sturm-Liouville problems proposed in 28 



In what follows we assume that / £ C 2 [— b, b], / ^ on [—&,&], /(0) = 1 and denote h := 
/'(0) £ C. Any such function is associated with an operator T^. For convenience, from now on 
we will write Tf instead of and the integral kernel of Tf will be denoted by K f . Together 
with the function / let us consider the function 1/f. Notice that if / is a solution of (2.30) 
with the potential qf = f"/f, the function 1/f is a solution of the equation 

+*//(*))« 

with the Darboux associated potential 

//V 

9i// = -?/ +2^ j J ■ 
Consider the transmutation operator T\jj which satisfies the equality 



= 



dx 2 



+ 9i// (a>) T 1/f [u] =T 1/f 



dx 2 



for any u £ C \—b, b] and transforms x into the functions ipk(x), k £ No defined by (2.24), i.e. 

..'<■_.,. i .. i r. .,. - r " (3 9) 



Ti//x = tpk{x) for any k £ N . 



It follows from (3.8), (3.9) and definitions (2.25), (2.26) and (2.39) that generalized wave 



polynomials are images of wave polynomials under the action of operators Tf and Ti/y, i.e. 

Tfp n = u n and Ti/fp n = v n for any n > 0. (3.10) 

In 26 explicit formulas were obtained for the kernel ~K\/f(x,t) in terms of I£f(x,t). In 
order to obtain a simpler expression for the integral kernel f{x,t) we assumed that the 
original integral kernel K/(x, t) is known in the larger domain than required by definition 
(3.2). Namely, suppose that the function JCf(x,t) is known and is continuously differentiable 
in the domain n : —6 < x < b, —b < t < b. It is worth mentioning that such continuation of 
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the integral kernel to the domain II is always possible due to (3.3), (3.4) and the general theory 



of Goursat problems. We refer the reader to 26 for further details. In Corollary 7.2 based on 



pseudoanalytic function theory we obtain another expression for the integral kernel Kx//(a;, t) 
suitable in the case when the integral kernel K/(a;, t) is known only in the natural domain 
1*1 < \x\ < b. 



Proposition 3.8 26 . Under the conditions of Theorem 3.6 the transmutation operators 
Tf and T\i f are related by the expressions 

T v/M(*) = ^ 



f(n)T s [u'\{ri)dn + u{0) 



and 



T f [u]{x) = f(x) 



o 7^ T i//["']^)*7 + w (°) 



valid for any u G C [— b,b], and their integral kernels Hf(x,t) := K/ (x, i; ft.) and K 1 /^(x, t) := 
JZi/f(x, t; —h) are related by the expressions 



K 1/f (x,t) 



1 



and 



K f (x,t) = -f(x) 



dt K f (s,t)f(s)ds + ^P-f(-t) 



1 f'fO) 



(3.11) 
(3.12) 



The following commutation relations immediately follow from Proposition 3.8 



Corollary 3.9 26 . The following operator equalities hold on C l [— b, b]: 



dwfTyt = fT f d x 
dx J Tf = J Tl/fdx ' 



Example 3.10 [26]. Take f(x) = x + 1 and some segment [—b,b] C (—1, 1). Then qj = 
and qi/f — (^jxp ■ For the potential qf the transmutation operator To is obviously the identity 
operator, hence K(a;, t; 0) = and from Proposition 3.3 we obtain 



K f (x,t) := K(x,t;l) = 1/2. 



Hence we get from (3.11 ) that 



K 1/f (x,t) 



t-1 
2(x + 1) ' 



It is easy to see that the integral kernels K/ and 'K.i/f satisfy the Goursat problems (3.5), (3.6) 
and relation (3.12) holds. 



For more examples of transmutation integral kernels we refer the reader to 26 . 



4. Goursat-to-Goursat transmutation operators 

Let u and u be solutions of the equations Ou — and (□ — q(x)) u — in R where R is 
a square with the vertices (±6, 0) and (0, ±6), respectively such that u = TfU. Consider the 
operator Tq mapping the Goursat data corresponding to u into the Goursat data corresponding 
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to u, 

u(x, x) \ / u(x, x) 

Jl(x,—x)J \u(x,—x) / 

Denote tp(x) := u(x, x), ip{x) := u(x, —x), := u(x, x), ty(x) :— u(x, —x). We have <p(0) 
■0(0) and <&(()) = ^'(0). The solution of the wave equation u is represented via its Goursat 

equalities 



Tn 



data as follows u(x,t) — (p( s ^) + — ¥>(0). Application of the operator T/ gives us the 



't + i N 
K/(x,tW — — )dr, 



T 'H £ i J )=*( £ i Li ) + >/(-.^( I i Li ) 



' x — t 



)dr 



and Tf[tp(Q)] = <p(0)f. Considering the value of these expressions when t — x and t = —x and 
introducing obvious changes of variables we obtain the following relations 



$(V) 



ip(x) 



and 



= p(0) 



Kf(x, 2i - x)y(t)dt + -0(0) 



K/(x, 2i + x)ip(t)dt + ip(x) 



K f (x, 2t + x)ip(t)dt - ip(0)f{x) 



K f (x, 2t - x)ip(t)dt - ip(0)f(x). 



Due to the boundedness of the operators Tf and TJ 1 as well as to the continuous dependence 
of the solu tion s of the Goursat problems under consideration on their respective Goursat data 
(see, e.g., 46, Sect. 15]) the operator Tq together with Tq 1 are bounded on vector functions 



from C 1 [—6, 6] x C 1 [— b, b] equipped with a suitable, for example, maximum, norm. 
It is easy to establish the following mapping properties of Tq, 

'u 2 n-i(x,x) N 
K u 2n -i(x,x) / 



(4.1) 



and 



(4.2) 



u 2n (x,x) 

< -U 2n (x,x) / 

Indeed, consider the wave polynomial p2 n -i(x,t), n — 1,2,..., see (2.39|. Its Goursat data 



have the form p 2n -i(x,x) = p2 n -i(x, — x) = 2 n ~ 1 x n . The image of p 2n -\ under the action of 
Tf is u 2n -i, see (3.10). Using the property ( 2.29| ) we obtain (4.1 1. Analogously, consideration 
of p2n leads to (4.2). 



Proposition 4.1. Let u be a regular solution of the equation (□ — q(x)) u = in R such 
that its Goursat data admit the following series expansions 



and 



^ oo 

- (u(x,x) +u(x,-x)) = c u (x,x) + ^ c n u 2n -i(x,x), 

n=l 

1 X 

- (u(x, x) - u(x, -x)) = ^ b n u 2n (x, x), 



both uniformly convergent on [—6, b]. Then for any (x, t) E R, 

oo 

u(x,t) = C Q U {x,t) + 2J {CnU2n-l(x,t) + b n U 2n {x 7 t)) 



(4.3) 



(4.4) 



(4.5) 
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and the series converges uniformly in R. 



Proof. Notice that if ([4^5]) is valid then (|43j) and hold. Indeed, due to fl2.29[ ) 

- (u(x, x) + u(x i —X)) = CqUq(x, x) 



^ oo 

+ ^ X! ( c ™ ( u 2«-i(a;,a;) + u 2n _i(x,:r)) + 6„ (u 2 „(a:, x) - u 2 „(x,x))) 

71=1 

oo 

= c o u (x,x) + ^ c n u 2n -i(x,x) 



and similarly (4.4) can be verified. 



Suppose that (4.3| and (4.4 1 are valid. Then we have 



u(x, x) 
u(x, —x) 



\ fu Q (x,x)\ ^ 

J =co U(, >a; )J + g 



U2 n -i(x,x) 
u 2n -i(x,x) 



From (J4TJ and ( |4.2[ ) we obtain 

=T -i( u(x,x) 
ip(x)J ' G \u(x,-x)J V I 



n=l 



u 2n (a;,:r) 
-u 2n (a;,a;) 



co + Er=i 2 " _1 (cn + &n)x" 

CO + E^ll 2 "" 1 ^-^)*" 



(4.6) 



Under the conditions of the proposition and due to the boundedness of Tq 1 we have that both 
series in the last equality are uniformly convergent on [—b,b]. Then applying Proposition 1 
from [18] we obtain that the unique solution u of the Goursat problem for the wave equation 
in R and with the Goursat data (^j) has the form 



u(x,t) = C p (x,t) + ^ ( c nP2n-l(x,t) + b n p 2n (x,t)) 



(4.7) 



and the series converges uniformly in R. Since u = TfU, Tfp n = u n and Tf is bounded we 
obtain the uniform convergence of the series (4.5 1. □ 



5. Transmutations for the Vekua operator 
By analogy with [5] we introduce the following pair of operators 

Vi=T f H+jT 1 /fZ (5-1) 

and 

V 2 =T 1/f K+jT f l. (5.2) 

Proposition 5.1. The following equalities hold for any M-valued, continuously differen- 
tiable function w defined on R. 

b\ - ^C^j V lW = V 2 (chw) , (d* + ^C^j V 2 w = V, (b\w) . 
d z - ^Cj Viw = V 2 (d z w) , (d z + ^-C^j V 2 w = Vi (d z w) . 
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Proof. The proof consists in a direct calculation with the aid of the commutation relations 
from Corollary |3.9| For example, for w = u + jv we have 

^-^f C ) ViW = Uf^(jK(V lW ) \ +jjd Y (fl(V 1 w)) 

= \ (fd w -jd t T f u+jjd x (fT 1/f v)-d t T 1/f v 

= \ (T 1/f d x u-jT f d t u+jT f d x v - T 1/f d t v) = V 2 (dgw) . 

□ 



An immediate corollary from Proposition |5.1| is the fact that the operators Vi and V2 



map the hyperbolic analytic functions satisfying (2.33) into the solutions of (2.4) and (2.10) 



respectively. Moreover, they map powers of the variable z into the corresponding formal powers. 



Proposition 5.2. For any z € R, n G No and a € B the following equalities are valid 
Vi [az n ] = Z {n) (a,0;z) and V 2 [az n ] = z[ n) (a,0;z). 



Proof. The proof consists in the observation that for a — a' + jb' and z = x + jt one has 



= (a'+j6')E 



•mj.77 
J 1 



and the result follows from the formulas ( |2.20[ ), ( |2.21 1, (2.22) by application of the mapping 
properties (f3T8|), flO). □ 



Notice that both operators Vi and V2 are bounded on the space of continuous functions. 
Indeed, consider 



|VH = * (|(Vn 
1 



(Vi«o 



= g (\K(v lW )+i(v lW )\ c + \n(v lW )-i(v lW )\ c ) 



T i/f v \ 



<\K(V 1 w)\ c + \l(V 1 w)\ c = \T f u\ l 

From the boundedness of the operators Tf and T\/f we have 

max |Vito| < M (max \u\ c + max \v\ c ) 

whereM = max{||T/|| , |[Ti/y || }. Since |u| c < \ (|w+| c + |w~| c ) and |u| c < \ 
we obtain max|Vitu| < 2Mmax|z«|. The boundedness of V 2 is proved analogously. 
The inverse operators Vj -1 and V^ 1 have the form 



w 



Vi X =Tj x 1l+3T-f f l 



and 



V^ 1 =T- 1 K+3T7 1 1 



26 



Theorem 



and clearly are bounded. For the explicit formulae for TJ 1 and we refer to 
10]. 

The following statement establishes a relation between the derivatives of hyperbolic analytic 
functions and the generalized derivatives of their images under the action of the transmutation 
operator. 



Proposition 5.3. Let w be a hyperbolic analytic function in R and W — Viu> be a 
corresponding solution of (2.4). Then whenever the corresponding derivative ofw exists there 
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exists a derivative of the same order in the sense of Bers of the function W and vice versa, and 
the following relations are valid 



= w [ - 



id V; 



(df n -^w) = W^ 2n ~ 1 \ n = l,2,. 



Proof. From Proposition |5.1| we have 

w = v 2 (d zW ) . 



(5.3) 



(5.4) 



W is a solution of the succeeding Vekua equation (2.10|. Denote W\ = W. Any solution of 



(2.10) is the image of a bicompl ex a nalytic function under the action of the operator V2, so 



W\ — V 2 W\. Due to Proposition [5Tj we have W\ = Vi (d z W\). Thus, W = Vi (d z w) because 
from (5.4 1 W\ — d z w. Now (|5.3|) can be easily proved by induction. □ 



6. Expansion theorem 



Theorem 6.1. A solution W of (|2.5|) can be represented as a uniformly convergent series 

(6.1) 



W(z) = u(z) +jv(z) = z(n) (on, 0; z) 



n=0 

in R iff the functions u(x,x) and u(x, —x) admit the uniformly convergent series expansions 
on [—6,6] of the form (4.3) and (4.4) where the coefficients of the expansions are related by 
the equalities 

a n =c n +jb„, ?i = 0,l,... 
and the coefficient 6q defines the value of v at the origin. 



b n := T(a n ), n = 0, 1, By (2.28) we have 



Proof. Suppose that W is a solution of ( 2.5 1 such that (6.1 ) holds. Denote c n :— lZ(a n ) and 



and 



Then 



Z^(l,0;z)=u 2n -i(x,t)+jv 2n (x,t), n= 1,2,... 
Z {n) {j,0;z) = u 2n (x,t) +jv 2n -i(x ) t), n= 1,2,... 



W(z) = c u (x,t) + b jv (x,t) + Y ( c « («2n-iOM) +jv 2n (x,t)) + b n (u 2n (x,t) + jv 2n -i(x,t))) ■ 

n=l 

Hence 

00 

U := K(W) = C U + Y ( c nU2n~l + 6„U 2 n) 
n=l 

and 

00 

v := T{W) = b a v Q + ^ ( c nV2n + 6n"2n-l) • 
n=l 

Thus, 

00 

u(x,x) = c u (x,x) + (c n u 2n -i(x,x) + b n U2n{x 7 x)) 



n=l 



20 



VLADISLAV V. KRAVCHENKO AND SERGII M. TORBA 



and 



u(x, -x) = c u (x,x) + ^2 (c n U2n-i(x,x) - b n u 2n (x,x)) . 



We obtain that (4.3) and (4.4 1 hold 



43 
4.5 



Now let us assume that W is a solution of (2.5 1 in R such that for u — 1Z(W) the equalities 
and (4.4 1 hold. Then due to Proposition 4.1 the function u admits a series expansion 



uniformly convergent in R and at the same time u — T ^ u is a solution of a Goursat 
problem for the wave equation with the Goursat data (^J defined by (4.6). u has the form 



Now, let us consider the Goursat problem (2.33), (2.35) with 



$(x) := 2(p(x) -CQ + bo and := 2ip(x) - c - b 

where bo € C is an arbitrary number. Thus, 

oo oo 

$(a:) = 2 " ( c ™ + M * n and *(*) = ^ 2" (c n - b n ) x n . 

n=0 n=0 

Due to Proposition 2.6 the unique solution of this Goursat problem has the form W(z) 
&nZ n where 

n=0 

a n = P + (c n + b n ) + P~ (c„ - b n ) = c„ + jfo„, n = 0, 1, 



Application of the operator Vi to W gives us a solution W of (2.5) in R in the form of a 
uniformly convergent series (6.1) with u — 1Z(W). □ 



Theorem 6.2. Let W be a solution of (2.5) admitting in R a uniformly convergent series 



expansion (6.1 ). Then there exist its derivatives in the sense of Bers of any order and the Taylor 



formula for the expansion coefficients holds, 



WM(o) 
n! 



(6.2) 



Proof. Under the condition of the theorem we obtain that w = *V 1 l W has the form ( 2.36 1 



with the expansion coefficients (2.38). The derivatives of w of any order exist because they 



reduce to the differentiation of power series expansions of $ and ^ from Proposition 2.6 



Application of the operator Vi to the even derivatives and of V2 to the odd, due to Proposition 



5.3 gives us the corresponding derivatives of W in the sense of Bers as well as the formula (6.2 ) 



□ 



Even if the conditions of Theorem 6.1 are not fulfilled and a solution of (2.5) cannot be 



represented globally as a uniformly convergent formal power series of the form (6.1), it may 



be arbitrarily closely approximated by finite linear combinations of the formal powers. The 
following analogue of the Runge approximation theorem from classical complex analysis is a 
corollary of the existence of the transmutation operator Vi . 



Proposition 6.3. Let W be a solution of (2.5) in R. Then there exists a sequence of 

J2n=o (°"> 0; z ) uniformly convergent to W in R. 



polynomials in formal powers P, 



Proof. Consider w = Vj 1 W. By Proposition 5.1 we have V2(dgw) — 0, hence w is a 



solution of (2.33) in R. Despite the domain R is smaller than the domain R, the proof of 
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Proposition |2.10] may be applied without changes to the function w. Hence for an arbitrary 
e > there exists a polynomial p n = X^Lq a n zn such that for any (x, i) G R 

£ 



\w(x,t) -p n (x,t)\ < 



iVill' 



where ||Vi|| < oo is the norm of the transmutation operator Vi. Applying Vi and Proposition 
15.21 we obtain that 



n=0 



< ||Vi|| • max_ \w(x, t) — p n (x, t)\ < e. 
(ij)eR 



□ 



7. Integral kernels as solutions of Vekua equations 

In Section[5]we showed that combinations (5.1 1 and (5.2) of the transmutation operators Tf 
and Ti/f are related with Vekua operators and formal powers. In this section we consider the 
combination 



K = K / -jK 1// 
of the integral kernels of the operators Tf and T\if . 



(7.1) 



(2.51. 



Theorem 7.1 . The function K given by ( 7.1 1 is a solution of the hyperbolic Vekua equation 



Proof. The proof immediately follows from relations (3.11 ) and (3.12) and conditions (2.6), 

□ 



(2.71 



The way of constructing the integral kernel T&Uf by the integral kernel K/ given by (3.11 ) 
requires the knowledge of the kernel K/(x, t) not only on a natural domain \t\ < \x\ < b 
sufficient for defining the transmutation operator, but on a larger domain \t\ < b, \x\ < b. 
Pseudoanalytic function theory allows us to obtain another method of reconstructing K-i/f by 
the known K/ and vice versa. Resulted formulae are simpler than those mentioned in 
Remark 21]. 
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Corollary 7.2. The integral kernels and ~K-i/f of the transmutation operators Tf and 
Ty i f are related by the expressions 



Ki//(x,£) = - 



1 



fWdtKfMdr, 



'* f'(x) 
o J\ x ) 



"* f'(0) 



(7.2) 



and 

K f (x,t) = -f(x) 



d t K 1/f (n,0)dr] 



\ K y f (x,Odt-^ 



'* f'(0) 
K 1/f (x,Odt+ J -±±f(x). (7.3) 

o L 
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Proof. Suppose first that qj :— € C 1 [— 6, b\. In such case we also have q 1 /f = ~ 
qj e C 1 [— b, b] and the integral kernels K/ and Kjyy are twice continuously differentiable with 



Hence we obtain from Theorems 



2.2 



and 



7.1 



respect to both variables and sati sfy e qua tion ( 3.5 ) with potentials qf and qi/f correspondingly 



for the function K t using the operator given 



by (2.14) that the function ~K-i/f differs from the function 



f( V )d t K f (r,,0)dr, 



* f'(x) 
dx K f (x,Odt+L^ 

o T\ x l 



by the term cf 1 , where c is a constant. From the condition Kx//(0, 0) = — ^p- 
Q) we find that c = 



(condition 



Suppose now that qj G C[—b, b]. We proceed as follows. Consider a sequence of polynomials 
{/n}neN such that /„(0) = 1 and /" -> f", f' n -> /' and /„ -t / uniformly on [-b, b] as n -t oo. 
See, e.g., [6j Proof of Theorem 11] for a possible construction. Since the function / does not 
vanish on [—6,6] we may additionally assume that all functions /„ do no vanish on [—6,6]. 
It is easy to see that q n :— j 1 - — > q uniformly and q n 6 C 1 [— b, b\. Every /„ is associated with 
the transmutation operators Tf n and 2\/y n with the integral kernels K^- re and K.i/f n . Consider 



functions H t (u, v) — K.f n (u + v, u — v; /^(0)). It is well known (e.g., 46 ) that the Goursat 



problem (3.3), (3.4) is equivalent to either of the systems of integral equations 



and 



(H(u,v) = | + Jo G(u',»)du', 
\G(U, V) = |ff(u) + f q{u + v')H(u, v') dv' 

f H( U , v) - | + \ H q(s) ds + $; F(u, «') d«', 
\F(«,«) = ^q(u' + v)U{u',v)du\ 

where G = d u H and F = ^H, from which we conclude that the functions Hy re , d u 1If n and 
d v Hf n are uniformly convergent to the functions Hy, d u llf and d v llf, respectively. Hence the 
partial derivatives of are uniformly convergent as well, and taking in ( |7.2| the limit as 
n — > oo we finish the proof. 

□ 



The second formula (7.3) easily follows from (7.2) by changing / to 1// 



Example 7.3. For the functions K.f and TZyf from Example 3.10 we get K = \ — j 2 ( x +i) 
and 



2 \ J 2(x + l) 2 2(x + 1) J 2{x + l) 2f(x) 



K. 



Formulae (7.2) and (7.3) can be verified as well 



8. Formulae for integral kernels of transmutation operators 



It follows from Theorems 6.2 and 7.1 and from (2.12) that the knowledge of the derivatives 
<9™K(0,0) allows one to obtain the generalized Taylor coefficients (6.2) for W = K. In this 
section we present several formulae for computing these derivatives whenever the derivatives 
q( m )(0), m < n — 1 are known. 

First we illustrate the expansion theorem with an example. 
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Example 8.1. For the function K from Example 7.3 by (2.12) we obtain 

1 



K[°](0,0) = i + ij, KM (0,0) 



I + Oj, K[" +1 l = jd t KM = 0, n > 1. 



(0,0) 



2(x + l) 

Hence the expansion coefhcients arc a,Q = | + ^j, ai = — g + Oj and a„ = for ?i > 2. 
The first two functions of the systems {vfc} and {i/'fe} are equal to 

1 



tpo = x + 1, 



"00 

0i 



2- + 1' 

a: 3 + 3x 2 + 3x 



3(x+l) ' 
the first two formal powers are given by 

{a, 0; z) = 71(a) ■ tp + jl(o) • Vo, 
Z« (a, 0; *) = ft(a) • + jt^o) + jI(o) • + jtpo) 
and we find that indeed, 



K(ar, t) = ZW (I + Ij, 0; z) + Z« (-i + 0j, 0; z) . 



We are able to calculate the expansion coefficients a n whenever we know values of the 
derivatives 9™K(0, 0). Recall that the function K/ is the integral kernel of a transmutation 
operator if and only if the function Hf(u,v) = K/(u + v, u — v) is the solution of the Goursat 
problem JO), dO). We have 



d?K f (x,t) = d?H f ( 



x + t x — t 



1 



(d u -d v ) n K f (u,v) 



x-{-t x — t 



Moreover, it follows from (3.4) that 



and 



d - Ufiu ' 0) -{y^\ u ), n>l, 




(8.1) 



(8.2) 



(8.3) 



Hence to calculate the value of 9™K/(0, 0) it is sufficient to transform the terms d^dy~ k Hf, 
where 1 < k < n — 1, into the terms involving only derivatives of the form dfjtlf or 9™H/ and 
some derivatives of the potential qj. Such transformation may be performed using equation 
(3.3) which allows us to reduce the orders k and n — k in the term d k d™~ k Hf by 1. After 
several applications of equation (3.3) and the Leibniz rule we obtain a finite sum of terms of 
the form 

q ( ; i) (u + v)..... q ^\u + v)di /v H f (u,v). 



(8.4) 



It is easy to see that the number I of factors q^ ni \u + v) satisfies the inequality 1 < l< 
min(fc, n — k) and that m + . . . + ri£ + 2£ + d = n. 
Moreover, consider a term of the form 



4) obtained as a result of application to H/ of 



some derivatives d u and d v , equation (3.3) and the Leibniz rule. If during this process we 
change every derivative d u by —d v and every derivative d v by — d u , we obtain the term of 
exactly the same form, with the only difference that the last derivative d^, H/ changes into 
d^ u Ilf, and the coefficient (—1)™ appears. Hence, we may group all the terms obtained after 
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the transformation of the expression (d u — d v ) n Ht(u, v) into the pairs of the form 

e 

(d d u + {-l) n d$)U f {u,v) ■l[q { p ) (u + v). (8.5) 

i=i 

Now we derive the recurrent relation for the derivative (d u — d v ) n Hf. We may parametrize 
each term of the form (8.5) by the following parameters. 



(1) A number £ of the factors q£^ {u + v). Since we need two derivatives to obtain one factor 
q(u + v) with the use of (3.3), the number £ satisfies 

n 



o<e< 



(8.6) 
(8.7) 



(2) A number d of derivatives of Hf satisfying 

< d < n - 21 

(3) An ordered sequence of non-negative numbers ri\ < ri2 < . . . < ni (orders of derivatives 
of the factors qf(u + v)) satisfying 

ni + ri2 + • ■ • + fi£ + 11 + d = n. 



(8.8) 



We denote by (— 1) S^-fni ne \ the coefficient at the term ( 8.5 ) in the final representation of 
(d u — d v ) n Hf after the described transformation. The factor (— 1)^ is used to make the number 
S? > / s non-negative. We assume that S? ,, , v = whenever either of the conditions 



.6)-(8.8) fails 



Lemma 8.2. The number of different lists of parameters £; d; (ni, . . . , ng) satisfying condi- 
for a fixed n coincides with the number p(n) of partitions ofn, i.e., the number 



-CI 



fcions 

of ways n may be represented as a sum of ordered positive integer terms, see, e.g., [lOf , [llf . 
The number p{n) is also known as the number of different Young or Ferrers diagrams. 

Proof. Consider a list of parameters £; d; (n\, . . . , ni) and put into correspondence to this 
list a partition 1, 1, . . . , 1, n\ + 2, n.2 + 2, . . . , ng + 2. Due to the condition (8.8 ) this partition 



is a partition of n. It is easy to see that the described correspondence is one-to-one. □ 

Note that (d u — d v )qr(u + v) = 0, hence in the case d = we have 
/ e 

(d u -d v )l (n f (u,v) + {-i) n n f (u,v)) n<zr> + «) 
^ i=i 

= (1 + (-1)") (d u + (~lT +1 d v )H f (u, v) ■ ft q^(u + v), (8.9) 

i=i 

and in the case d ^ we have 

(d u - d v ) ( (d d u + (-l) n b^)Hf(u, v) n q { p\u + v) 
v i=i > 

= + {-i) n+i d d v +i )n f {u,v) - (dt 1 + (-i) n+i ar i )( (Z/ («+t')H / K«))) n^f^") 

■i=i 

= (d d u +1 + (-ir +1 d d v +1 )H f (u, v) ■ n qf i] (u + v) 

i=i 
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fc=0 



?? } (« + «)(^ 1 - fc + (-l) n+1 ^- 1 - fe )H / (u,t;) • n?/" ^ + v )- ( 8 - 10 ) 



Also note that the coefficient at the term (d™ + (— l) n d™)Hf(u, v) is equal to 1/2 for n — and 
1 for n > 1, hence S° . () = 1/2 and = 1 for n > 1. From the relations (JsTsJ) and ( |8.10 | 

and the latter initial condition we obtain the following statement. 



Proposition 8.3. The coefficients <S?t) ^ 



with the hsts of parameters t\ d\ (n\ , . . . , ng ) 



satisfying conditions (8.6 1-(8.8 1 can be calculated by means of the following recurrent relations 

S l-X(ni,...,n e ) ^l 1 + i^ l T) S i:d~l;(n 1 ,...,n e ) 

fd + n k \ (8.11) 

/ , I n fe / •«- 1 ;d+n*+li(ni,...,nfc-i,n*+i,...,n«) " » — i, 



on+l nn 

D e-,d;(n 1 ,...,n e ) ~ °£;d-l;(n 1 ,. . .,n e ) 

+ E 

n k e{n ly ...,n e } 



d + n k \ (8.12) 



where _ r i means that only different values of n,- are used in the sum, with the initial 

conditions 



(8.13) 



The following proposition presents a direct formula for the coefficients S™. d ., ^ 
involving recurrent relations. 



not 



PROPOSITION 8.4. The following formula holds for the coefficient S™. d .^ ni 
and the list of parameters I; d; (m, . . . , rig) satisfies conditions (8.6 )-(8.8 ): 



where £ > 1 



d di+o-i + l d f _i+(Tf_i + l 

S e ; d;(ni,...,n £ ) = E EE'" E 

(<ri,...,<Ti)e©(ni,...,n*) di=0 ^2=0 di=0 



H (* + ^ (1 + (-ir+-+-^)(l - 5(d)5(i - 1))), (8.14) 
i=l v CTi / 

where (oi, . . . , at) € &(ni, . . . , ng) means that the sum is taken over all distinct permutations 
of the numbers m, . . . , nt, and 8{t) = 1 if t = and S(t) — otherwise. 



Proof. To arrive from the term S] 



to a term 5q.-.q using recurrent relations 
111) and (|8.12[) we have to perform exactly € times the step of removing one of the n^'s from 



£;rf;(ni , . .. ,n^) 



the list (m, • ■ • j H^). Denote by dj, . . . , the values of the parameter d on each of these steps 
and by cti, . . . , ct£ the number rik removed from the list (ni, . . . , rig) on the corresponding step 
1, . . . , I. Then the numbers d\, . . . , dg satisfy the following inequalities 

< di < d, 

<d 2 < 1 + di +cri, 



< di < 1 + c?£- 1 + <x^_ i ■ 
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Each distinct permutation of the numbers n\, . . . ,ri£ and each list of parameters d\, . . . , di 
give us a different way of getting from the term S2. d .r ni n( \ to the term Sq.~.q. Note that 



the factor 1 + ( — 1)™ for some n' < n appears in (8.11) only in the case when we cross the 
coefficient S^,. .^ n , n , ) having d' = 0. Note also that if d = initially, then d\ — and only 



12 ) is applicable, hence the factor 1 + ( — 1)" does not appear in this case. So the factor 
1 + (—1)™ appears only in one of the following cases: 1) d ^ and d\ — 0, 2) di ~ 0, i > 2. 



It follows from condition (8.8 1 that in each of these cases n' = 2£' + n[ + . . . + n' t , therefore 
(— 1)™ = (— +«£'. Hence we conclude that the following relation holds 

Sl M nu..,n t) = E E ( dl ^ ai )(l + (-ir + ' + --W(l-^))) 

(o-i,...,CTf)£e(rai,...,7v) di=0 ^ 1 ' 

x E ( d2+(T2 )(i + (-ir + - + ^(rf 2 ))--- E ( de+<7e )(i + (-irs( de) ), 

d 2 =0 V CT2 / ^=0 \ * / 

where the first sum is taken over all distinct permutations of numbers n\,...,ni, and the 
function S(-) is used to the encode mentioned above conditions 1) and 2). Finally note that the 
conditions 1) and 2) may be jointly represented by the expression S(di)(l — 5(d)8(i — 1)). □ 

As a result we obtain the following representation for the derivative (d u — d v ) n Hf(u, v), 
n > 1. 



Proposition 8.5. Let n > 1 and q f e C^-^-fr, 6]. Then 



- o ) n H/(u,«) = + (-ird:)n f ( u ,v) 

[n/2] n -1l I 

+ E^ 1 )' E( 5 « + (- 1 )" 9 ") H /( W ^) E ^(n I ,..,n,)ri5/ n<) ( U + U ) 



i=0 



ni + ...+n^— n— 2^— 
0<ni <. . .<n^ 



where the coefficients Sf. d ., ni % are given by Proposition 8.3 or Proposition 



depend on qf. 



and do not 



The following proposition is a corollary of (8.1|, (8.2), (8.3) and Proposition 8.5 



Proposition 8.6. Let n > 1 and q f G C( n_1 > [-&,&]. Then 

[n/2] n -2l 

2 n 



-, / [n/2\ n-2l 



d=0 



e ^;(^...^)n^ i) (°))' ( 8 - 15 ) 



ni + ...+n^— n— d— 2^, 
0<ni<...<ne 



where we set q jTj (0) := = /'(0), the coefficients S". d ., ni n s are given by Proposition 
or Proposition 8.4 and do not depend on qf, and 5(d) = 1 if d = and 6(d) = otherwise. 
Formula (8.15 1 aJso holds for n = for any e C[— 6, 6]. 



8.3 



As a consequence of Theorems |6. 1[ [6~2] and [7T] we obtain the following representation of the 
integral kernel K/. 



CONSTRUCTION OF TRANSMUTATION OPERATORS 



27 



Theorem 8.7. Let qt g C°°[— b, b] be a complex-valued function and f be a particular 
solution of (2.30) such that f ^ on [— b, b] and normalized as /(0) = 1. Denote h := 
/'(0) € C. Suppose that the functions g\{x) := ^(K/(.x, x) + Kf(x, —x)) = 5 + |J ?/( s ) ds 
and g 2 (x) := |(Kj(i,i:) — K^(x, — 2;)) = | <7j(s) <is admit uniformly convergent on [—0,6] 
series expansions 

00 

9l(x) = C U (x,x) + ^ CnU 2n -l(x,x) 



and 



92{x) = ^2 b n U 2 n(x,x). 



Then the coefficients {c n }„>o and {b n } n >i may be found by the formulas 

W) . d?K 1/f (0,0) . 
c n = —. , o n = —. , if n is even 



and 



a t "Ki /f (0,0) ^"K/(0,0) ., . 

. , b n = J - , if n is odd, 

n\ n! 



where the derivatives <9™K/(0,0) and 9"Kx//(0,0) are given by (8.15) for the potentials qj 
and qi/f = —qf + 2(/'//) , respectively. Also, for any (x,t) £ R 



K f (x,t) = c Q u Q (x,t) + ^ {c n u 2 n-i(x,t) + b n u 2n (x,t)) , 



(8.16) 



and the series converges uniformly in R. 



Example 8.8. For the potentials qf and q^/f from Example 3.10 we have 
q f - 1 \0):=f(0) = l ) q f n) (0) = 0, n>0, 



and 



g y ) (0):=(l//)'(0) = -l, g™(0) = 2(-l)> + l)!, n > 0. 

the first coefficients A , ^ for n < 6 have the following values 

^0;l;() = 1' 



By Proposition 



8.3 



c2 _ i o2 
°0;2;() — x i °1;0;(0) 



^0;3;() — 1) ^1:0:(1) — 1> ^1;1;(0) ~~ ^' 



c4 

D 0;4;() — 
c5 _ -I . o5 
D 0;5;() — X > J l;0;(3) 



1 l ? 5 



0;(2) 

(2) 



<? 4 

J l;l;(l) 



C,4 

°1;2;(0) 



5. 



2;0;(0,0) 



^' ^1;2; 



s, 



0;6;() 



= l; 



sr. 



1;0;(4) 



= 1, 



c-6 

°2;0;(0,2) 



10, 



06 

J 2;0;(l,l) 



06 

D i;i;(3) 
5, 



(i) 
= 4: 



^' ^l;3;(0) 



^1 ^2:0;(0,1) 



^' ^2;1;(0,0) — 



06 

i;2;(2) 



- H) ^1;3;(1) — 9 > ^1 



c6 



(0,1) 



We can check that formula (8.15) works, e.g., 
1 

64 



15, i$2;2;(0,0) — ^3;0;(0,0,0) 



24) 



4;(Q) — ^' 
10. 



5 t 5 K l7/ (0, 0) = i (240 + (-1) 1 (1 • (1 + (-1) 5 ) • 48 + 5 • 24 + 5 • 16 + 5 • 

+ (-1) 2 (6 • (1 + (-1) 5 ) • 8 + 10 ■ 8)) = 0, 
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and 



a?K v/ (Q,0) 



^8 (-1440 + (-1) X (1 • (1 + (-1) 6 ) • (-240) + 4 • (-96) + 11 • (-48) 



+ 9 • (-48) + 6 • (-96)) + (-1) 2 (10 • (1 + (-1) 6 ) • (-24) + 5 • (1 + (-1) 6 ) • (-16) 

+ 15 ■ (-16) + 15 ■ (-16)) + (-1) 3 (10 ■ (1 + (-1) 6 ) ■ (-8))) = 0. 



Now consider an example when all the coefficients b n , c n from (8.16) may be easily calculated 



explicitly from (8.14) and (8.15) 



Example 8.9. Consider qf = c 2 for some constant c taking the function / = e cx as a non- 
vanishing solution of the equation /" — <?// = 0. Then the Darboux transformed potential has 
theform^ =2(/'//) 2 -< 7/ = c 2 . 

Consider the function H/ from (8.1). Equation (3.3) for qf — c 2 reads as d u d v Tlf = c 2 Hf, 
hence 



(d u d v ) n H f (u, v) = £(-l) n - fe (? K^ n ~ fe H/K<;) 



k=0 

(d:+(-ird:)n f (u,v) 

[n/2]-l 

/■ ..'2/. / " \ g 



£ (-i) n ~ fe c : 



k=l 
n-1 



i-2k 



Kf(u,v) 



(8.17) 



+ £ (-i)"" 

fc=[n/2] 



k c 2(n-k) 



df- n K f {u,v). 



It follows from (8.2) and (8.3) that at the point u = v = every term in (8.17) involving a 
derivative of H/ of order at least 2 with respect to u or at least 1 with respect to v equals zero. 



Hence from (8.1 1 for any even n = 2m we obtain 

1 



9^(0,0) = ^(-irc 



. 9m\ f_l\m 2m+l /n 



in 



and for any odd n = 2m + 1 we obtain 



1 



22m+t 



in 



dr +i K f (o,o) = ^{-iTH 2m + 1 )a.H / (o,o) = ^ mc ^ J <1»> + 1 



22m+2 



(8.18) 



(8.19) 



Similarly for the kernel Kjyj we have 

(-l) m c 2m+1 (2m 



d 2m K 1/f (0,0) = -- 



c) 



2m +1 



2 2m+l y m 

Hence by Theorem |8 . 7| the expansion coefficients have the form 



Ki//(0,0) - 22m+2 



C2rn+1 



(-l)" 1 ^™^^ 
2 2m+l( m! )2 ' 

^_iyn+l c 2m+2 

2 2m + 2 m\[m + l)V 



(-l)" 1 ^™^^ 
2 2m+l( m !)2 ' 

(_l) m c 2m+2 

2 2m + 2 m!(m+ 1)! 



Now we show that it is possible to obtain the same coefficients from (8.14) and (8.15). Note 



that any term in (8.15) involving a positive order derivative of qf is equal to zero. Hence 
n% = n.2 = • • • = nt = and n = 2£ + d. Moreover, c? < 1 and is equal either to or 1 depending 
on the parity of n. So we have to calculate only one coefficient, S^ o ., q n for n = 2m and 
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S 1 !™ m for n = 2m + 1. Formula (|8. 14b for the coefficient S* 2 ™.^ m takes the form 



i;l;(0,...,0) 



J m;0;(0,...,0) 



°m;0;(0,...,0) 



1 d 2 + l 

d 2 =0 d 3 =0 d m =0 



We claim that 



m,-fc + l d m -l+l 

£ (l + *(dm-fc+l))-" E (* + *(*»)) = 
-fe+i=0 d m =0 



2fc + 1 + cZ m _ fc 
fc 



.20) 



Indeed, for = 1 we have =0 + + ^(^m)) = ^m-i + 3 = ( 3+d ™ x ) . Assume that the claim 
is correct for some fc, then for k + 1 we obtain 



-is-i+i 



2/c + 1 + d m _ fc 
k 



£ C 1 + <5(<*m-fc)) 

d m _ fe =0 

2A: + 1^ + /2fc + 2^ + + /2fc + 2 + d m _ fc _i s 
combinatorial identity = (») + ("7 1 ) + . . . + ("7 s ) + see, e.g., [TT 



and the induction step follows from the equality 2 ^ k j , y k ) — \ k+1 



we obtain from ( pT20| ) that S^™ 



For fc = rn — 1 



;(o,...,o) 



f-" 1 ) and d 2 ™ K/ (0,0) = ^cf™" 1 ) (c 2 ) 



T /2m\ 

2 2 m +l (, m J 



Similarly, for n = 2m + 1 we obtain from ( |8.20[ ) for k = m that S^t 1 



a t 2m+1 K/(0,0) 



- MO!" J2/2m+l\ / 2\m _ H 



(o,...,o) = ( 2 ™^) and 



2 2m+2 ^ m + 



2 2m + 2 



( 2 ^ +1 ), exactly as in ([0§ and J87l9] >. 



In the following example wc calculate the expansion coefficients b n and c n , n < 22, for 
a reflectionless pote ntia l in t he one-dimensional quantum scattering theory, q SO ch(x) = 1 — 
2sech 2 x, see, e.g., 26 , 30 . For this potential to the difference with from the previous 
examples neither the higher derivatives of qf nor the higher derivatives with respect to t 
of the transmutation kernel do not vanish identically. 



Example 8.10. The potential q SO ch(x) = 1 — 2sech 2 x can be obtained as a result of the 
Darboux transformation of the equation u" = u with the potential g C osh = 1 with respect to 
the solution f(x) = cosh a;. The transmutation operator for the operator A cos h = d 2 — 1 was 
calculated in [61 Example 3]. Its kernel is given by the expression 



K cosh (iE,t) 



1 Vx 2 -t 2 h{Vx 2 -t 2 ) 



x — t 



where l\ is the modified Bcssel function of the first kind. Even though in 26 a formula for the 



integral kernel K scc h was presented, it is not well suited for calculating higher order derivatives 



with respect to t. Using (7.2 1 we obtain another representation with the help of Maple 12 
software 



K scc h(a;, t) = — ( Ii(x) — Iq(x) tanhx + tanhrr 



Vx 2 



J -hWx 2 -s 2 ) 



ds 



'* (xs - x 2 )I (Vx 2 - s 2 ) + \Jx 2 - s 2 Ii(Vx 2 — s 2 ) 



ds 



.21) 



o (x-s) 2 

Despite the integral kernel K scc h could not be evaluated explicitly, it is possible to expand it into 
a Taylor series at the point t = 0. Expanding the kernels K cos h and K sec h into corresponding 
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Taylor series and taking a limit as x — > we obtain (again with the help of Maple 12 software) 
that 

K cosh (0,i) = 1 1- — t 3 + — t 5 - t 7 + 1 i 9 1 t 11 

coshy , i 4 32 76g 3g864 2949120 353894400 

+ t 13 - t 15 + - t 17 

59454259200 13317754060800 3835513169510400 

~ 1380784741023744000 t + 607545286050447360000 1 + °^ ' 

and 

K (0 t) = - - t + — t 3 - 1 < 5 + 1 {! 1 +9 j 1 xll 

sc v ' ' 4 96 3840 258048 26542080 3892838400 

1 . t 13 + —J—— t 15 - - — — i- — — t 17 



772905369600 199766310912000 65203723881676800 

' i 19 1 i 21 +o(t- 22 ) 



26234910079451136000 12758451007059394560000 
According to Theorem |8.7| we obtain the following table of the expansion coefficients for the 
integral kernel K scc h. Since all even coefficients are equal to 0, we do not include them in the 
table. To save the space, we represent all odd coefficients in the form of a fraction 2T1 +^ , . 



n 


l 


3 


5 


7 


9 


11 


bn 


i 

2 2 1! 


i 

2 4 3! 


2 
265! 


5 
2»7! 


14 
2 1 <>9! 


42 
2 12 11! 


(' n 


1 

2 2 1! 


3 
2 4 3! 


10 
2»5! 


35 
2»7! 


126 
2^9! 


462 
2 12 11! 



n 


13 


15 


17 


19 


21 


b n 


132 


429 


1430 


4862 


16796 


2 14 13! 


2 16 15! 


2 18 17! 


2 20 19! 


2 22 21! 


('■11 


1716 


6435 


24310 


92378 


352716 


2 14 13! 


2 16 15! 


2 18 17! 


2 20 19! 


2 22 21! 



For / = cosh a; we have h = f'(0) = 0. The derivatives at the point x — of the potential 
qi = 1 are equal to gj O) (0) = 1, ^ n) (0) = 0, n > 1, the non-zero values of the derivatives of the 
potential q(x) = 1 — 2sech 2 (a;) at the same point are given in the following table. 



n 


2 


4 6 


8 10 


12 


gW(Q) 


-1 4 


-32 544 


-15872 707584 


-44736512 




n 


14 


16 


18 


20 


q^(0) 


3807514624 


-419730685952 


58177770225664 


-9902996106248192 



The table below presents the values of the expression in parentheses appearing in 
(8.15), obtained from the derivatives of the potentials presented above with the coefficients 
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^7-d-(ni n t ) calculated by Proposition 8.3 All calculations were performed with the help 
of Matlab 7 and the results agree witrTThe ones obtained from the explicit formulae for 
the integral kernels (see the table for b n and c„ above). The value n = 21 is chosen only 
to keep the size of the table reasonable. Our Matlab program calculated exact values of the 
derivatives 9™K cos h(0, 0) for n < 30 and exact values of the derivatives <9"K SGC h(0, 0) for n < 26. 
To calculate the derivatives for larger values of n one have to use arbitrary precision arithmetic 
to overcome the limitation of the machine precision. 



n 


l 


3 


5 


7 


9 


11 


2 n+1 d?K cosh 


l 


-3 


10 


-35 


126 


-462 


on+1 nniz 


-l 


1 


-2 


5 


-14 


42 





13 


15 


17 


19 


21 


2" +1 a t n K cosh 


1716 


-6435 


24310 


-92378 


352716 




-132 


429 


-1430 


4862 


-16796 



9. Approximate construction of integral kernels 

In this section we discuss two methods of approximate construction of integral kernels of 
transmutation operators. The first method is based on the expansion theorem and computation 



of the generalized Taylor coefficients by Theorem 8.7 and formula (8.151. The second method 



is based on the completeness result of hyperbolic formal powers and of generalized wave 



polynomials (Proposition 6.3 and [18| Theorem 22]) and on approximating the Goursat data for 
the integral kernel by generalized wave polynomials. However we present only the description 
of the main ideas of these approximate methods and do not perform detailed analysis. 
Such analysis, involving studying of approximative properties of the functions {u n (x,x)}, 
convergency rate estimates and detailed comparison with existing numerical techniques, e.g., 
successive approximation and the series expansions from [3], goes beyond the scope of the 
present article. The authors hope to perform such analysis and present the results shortly. 

The first method is applicable in the case when we know derivatives of the function / G 
C n [— b, b) at the point x — of all orders up to some order n, n>2, i.e. when we know 
coefficients in the Taylor formula 

n 

f(x) = J2f^ k +o(x k ). 

fe=0 

Suppose additionally that the function / does not vanish on [—b, b] and normalized as f(0) = 1, 
i.e. /o = I- Then we can calculate Taylor coefficients up to the order n of the function 1/ /. For 
example, as it follows from |10[ Theorem 11.7] or directly from the Faa di Bruno formula |10| 
Theorem 11.4], 

1 ™ 

J v ; fc=i 



32 



VLADISLAV V. KRAVCHENKO AND SERGII M. TORBA 



where 



fk 



E 



mi + ...+m fe 



mi +2m2 + . . .+fcrrifc —k 
mi>0,... ,m& >0 



mi, . . 



+ mjt 
,m fc 



li/; 



(9.1) 



and (™i+-+"M = are multinomial coefficients. 

Having coefficients and fk, k < n, we easily calculate the Taylor coefficients up to the 
order n — 2 of the potentials qj = f"/f and q-yij = f ■ (1//)". By (8.15) we obtain the values 

we find the 



8.7 



of the derivatives 9J"K/(0,0) and 9J n K 1 /y(0, 0), m < n — 1, and by Theorem 
coefficients and c^, k < n — 1 in the representation ( 8.16[ ). The obtained truncated series is 
an approximation of the integral kernel Ky. 

If on the opposite the Taylor coefficients of the potential q/ and the value of the parameter 
h = f (0) are known, the reconstruction of the Taylor coefficients of the solution / also presents 
no difficulty. 



Example 9.1. Consider the function / = coshx as in Example |8. 10 One can easily verify 
that the first Taylor coefficients for the function 1/f calculated by (9.1) coincide with the 
coefficients in the known expansion 



1 



cosh: 



= sech x = 1 — 



1 



24 



61 
720 ^ 



+ 277 x 8 + Q(x 10 ) 
+ 8064 X +U[X h 



and the calculated Taylor coefficients for the potentials qf and q±/f coincide with the values 
presented in Example |8.10| 




We compare the exact integral kernel K scc h given by (8.21) with its approximations based 
on the generalized Taylor coefficients calculated in Example 8.10| To illustrate the dependence 
of the approximation error on the domain R and on the number of coefficients used, we present 
the resulting approximation error for 3 different domains corresponding to values 6 = 1,6 = 2 
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and 6 = 4, and for different number of generalized Taylor coefficients. The functions ipk, defined 
by ( 2.23 1 and necessary for constructing the generalized wave polynomials, are calculated using 
two Matlab routines from the Spline Toolbox: on each step the integrand is approximated by 
a spline with 5000 knots using the command spapi and then it is integrated using fnint. The 
resulting kernel K scc h and its approximations were computed on the mesh of f 00 x f 00 equally 
spaced points in R. Obtained results together with the graphs of the integral kernel K sec h and 
the approximation error are given below. 



6 = 1 6=2 6 = 4 



N Approx. error N Approx. error N Approx. error 



f 


0.12833 


5 


0.21204 


13 


1.8261 


3 


0.021458 


7 


0.043909 


15 


0.39987 


5 


0.0017866 


9 


0.0059553 


17 


0.070023 


7 


8.9155 • 10~ 5 


11 


0.00057228 


19 


0.010037 


9 


2.9655 • 10~ 6 


13 


4.1076 • 10~ 5 


21 


0.0011998 


11 


7.0469 • 10~ 8 


15 


2.288 • 10~ 6 


23 


0.00012146 


13 


1.2562 • 10~ 9 


17 


1.0182 • 10~ 7 


25 


1.055 • 10~ 5 


15 


7.3683 • 10~ n 


19 


3.7047- 10- 9 


27 


7.9493 • 10~ 7 


17 


7.3991 • 10~ n 


21 


3.3373 • 10~ 10 


29 


5.2453 • 10~ 8 


19 


7.3989 • 10~ n 


23 


3.3373 • 10~ 10 


31 


3.0562 • 10~ 9 




Figure 2. The approximation error of the integral kernel i4T S0C h in the domain R given by 
b = 4 by the generalized wave polynomial of order N = 31. 
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In order to estimate the error of approximation the exact formula (8.21 1 for the kernel K 



^-scch 



was used, where the corresponding integrals were computed numerically. As can be seen from 



.21 ) the integrands for t close to x involve subtraction of close values and division by values 



close to zero. Due to this fact the kernel K scc h computed in this way by (8.211 is not exact. 
This explains why in the above tables after a certain value of N the reported approximation 
error does not decrease, it is limited by the computational error of K soc h. 

To the difference of K sec h the integral kernel K cos h does not involve integration and may be 
computed more precisely. The approximations of the integral kernel K cos h were computed for 
b — 2 and the resulting absolute errors obtained for different numbers of expansion coefficients 
are presented in the table below. 



6 = 2 



N 


Approx. error 


N 


Approx. error 


N 


Approx. error 


1 


1.7878 


11 


0.00074903 


21 


1.316 • 10~ 10 


3 


1.088 


13 


5.2024- 10~ 5 


23 


3.3386 • 10~ 12 


5 


0.33724 


15 


2.8243 • 10~ 6 


25 


7.5051 • 10~ 14 


7 


0.063779 


17 


1.2312 ■ 10~ 7 


27 


6.9944 • 10~ 15 


9 


0.0081416 


19 


4.4042 ■ 10~ 9 


29 


6.6613 • 10~ 15 



The first method of approximate construction of the integral kernels requires the knowledge 
of higher order derivatives of the potentials q f and q\j f in the origin. In the case when 
such derivatives are unknown or hard to obtain analytically, instead of performing numerical 
differentiation we can approximate the integral kernel by approximating its Goursat data given 
by (3.6) with truncated series of the form (4.3 1 and (4.4). 

Indeed, suppose that the numbers Co, . . . , cjv an d 6i, . . . , 6jv are such that 

A? 



-(Kf(x,x) + K f (x,-x)) - c u (x,x) - ^ c n u 2n -i{x,: 



< £l 



and 



1 N 

- (K/(ar, x) - Kf(x, -x)) - ^ ^^(i, x) 



< £2 



for every x £ [~b, b]. Consider the functions 



N 



K(x,t) = c u (x,t) + 2J c n u 2n ~i(x,t) + ^ b n U2n(x,t) 



(9.2) 



and K/ = T* X K/ and K = Tj 1 K. Then by the definition of the Goursat-to- Goursat 
transmutation operator 

K f (x,x) \ T i / K f (x,x) \ 

and ( ^( a: ' a; ) ) = T^ 1 ( ^{x,x) 



K f (x,-x) G \K f (x,-x) 



hence due to the boundedness of the operator T r 1 



K(x,-x) G \K{x-x) 



max! max \Kf(x,x)—K(x,x)\, max |Kj(a;, — x) — K {x, — x) 



K x€[-b,b] 



xE[-b,b]' 
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< 1 1 1 1 1 max max K/(a;, x) — K(x, x)\, max |Ky(:r, — x) — K{x 1 — x)\ 

b,b] x£[— b,b] 

\ (K/(x, x) + K f {x, -x)) - l - (K{x, x) + K(x, -x)) 



< IIT^ 1 !! max 
~ 11 G " xe[-b,b] 



\ ( K f(x, x) - K f (x, -x)) - i (K(x, x) - K(x, -x)) 



< ||T G 1 ||( £l+£2 ), 



where we have used equalities | (K(x, x) + K(x, —x)j = cquq(x, x 
\ (K(x, x) — K(x, —a;)) = Yl n =i bnU2n(x, x). We obtain from the proof of 
for every £ R 

-i-i 



18 



x) and 
Theorem 31 that 



K^ar.t)-^,*) <3||Tc i ||(£i+£2), 



hence for every (x,t) e R 



\K f (x,t) - K(x,t)\ < 3\\T f \\ ■ HTg 1 !!^ +e 2 ) 



That is, if we approximate the function gi(x) := |(Ky(x, x) +K.f(x, —x)) = | + | Jg Qf( s ) ds 
by the functions {uo(x, x)} U {u2 n -i(x, x)} n >i and the function g2{x) := ^ (K.f(x,x) — 
K.f(x,—x)) = \ L gf (s) ds by the functions {u2 n (x, x)} n >i, we obtain an approximation 
having the form ( |9.2[ ) of the integral kernel. It is worth to mention that we do not need to 
know the Darboux associated potential qi/f for this method and we do not impose additional 
smoothness conditions on the potential qj. The approximation of the functions g\ and g 2 by the 
corresponding combinations of the functions u n (x, x) can be done in several ways. We may use 
the least squares method to obtain a good, however far from the best, approximation. We may 
reformulate the approximation problem as a linear programming problem and solve it. Under 
the additional assumption that the corresponding systems of the functions u n (x, x) satisfy the 



Haar condition, we may use the Remez algorithm. We refer the reader to 18 Section 6] and 
to (35j §7], (39} Chapter 6] and (40j Chapter 1] for further details. 



Example 9.2. We approximated the integral kernels K soc h and K cos h from Example |9.1 
by the second described method. All the parameters are taken as in the previous example. 
The Goursat data were approximated by corresponding functions u n (x 1 x) with the help of the 
Remez simple exchange algorithm, see e.g., 35 §7]. 



b = 2, K soch kernel 



N Approx. error 



5 0.0045993 

9 9.3687 -10- 6 

13 4.1549 -10- 9 

17 3.3354 • 10~ 1C 



6 = 4, K scch kernel 



N Approx. error 



13 7.4042 • 10~ 5 

17 2.342 • 10~ 7 

21 2.789 • 10~ 10 

25 4.9467 • 10~ n 



6 = 2, K cosh kernel 



N Approx. error 



5 0.0052907 

9 1.2563 -10- 5 

13 6.6227 -10- 9 

17 1.1813 -10" 12 

19 1.0325 • 10-" 



References 

1. H. Begehr and R. Gilbert, Transformations, transmutations and kernel [unctions, vol. 1-2 (Longman 
Scientific & Technical, Harlow, 1992). 

2. L. Bers, Theory of pseudo-analytic functions (New York University, 1952). 

3. A. Boumenir, 'The approximation of the transmutation kernel', J. Math. Phys. 47 (2006) 013505. 



36 



VLADISLAV V. KRAVCHENKO AND SERGII M. TORBA 




-2 -2 
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is obtained with the use of the Remez simple exchange algorithm. 
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